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Abstract 



We develop the perturbative QCD formalism for inclusive semileptonic B meson 
decays, which includes Sudakov suppression from the resummation of large radia- 
tive corrections near the high end of charged lepton energy. Transverse degrees of 
freedom of partons are introduced to facilitate the factorization of B meson decays. 
Ambiguities appearing in the quark-level analysis are then avoided. A universal dis- 
tribution function, arising from the nonperturbative Fermi motion of the b quark, is 
constructed according to the heavy quark effective field theory based operator prod- 
uct expansion, through which the mean and the width of the distribution function 
are related to hadronic matrix elements of local operators. Charged lepton spectra 
of the B — > X u tv decay are presented. We find 50% suppression near the end point 
of the spectrum. The overall suppression on the total decay rate is 8% for the free 
quark model, and is less than 7% for the use of smooth distribution functions. With 
our predictions, it is then possible to extract the Cabibbo-Kobayashi-Maskawa ma- 
trix element \V u b\ from experimental data. We also discuss possible implications of 
our analysis when confronted with the rather small observed semileptonic branching 
ratio in B meson decays. 



1 Introduction 



The studies of semileptonic decays in heavy mesons within the framework of per- 
turbative quantum chromodynamics (PQCD) dated back to the seventies 0. The 
mass Mq of a heavy quark Q provides a rationale of this approach. The advan- 
tage of the PQCD formalism is that it provides a natural normalization of decay 
amplitudes. This point is of great importance to obtain a model independent ex- 
traction of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, which are the 
key phenomenological parameters in understanding the symmetry breaking physics 
of the standard model. We have shown that PQCD is applicable to the exclusive 
B — > niv and B — > piu decays [[J], and an upper limit of \V u b\ around 2.7 x was 
extracted from experimental data directly. Recently, there have been controver- 
sies concerning the inclusive semileptonic branching ratio of B meson decays @, ||. 
Even new physics was proposed 0] to resolve the discrepancy between theoretical 
predictions and experimental outcomes. It has been argued that PQCD radiative 
corrections may play an important role in the semileptonic decays B — > Xtv |J. In 
view of this, it is essential to first sort out the correct PQCD contributions to these 
decays. 

In the pioneering works of Chay et al. |J and Shifman et al. M a systematic 
expansion of relevant hadronic matrix elements in semileptonic decays in the in- 
verse power of Mq was obtained by combining heavy quark effective field theory 
(HQEFT) and the method of operator product expansion (OPE). It was shown 
that the differential decay rate dT / dq 2 dEi can be consistently calculated, when it 
is suitably averaged over the charged lepton energy Eg. Here q 2 = (p£ + p u ) 2 is 
the lepton pair invariant mass with pi and p v the charged lepton and neutrino mo- 
menta, respectively. At leading order, the expansion reproduces the naive parton 
model predictions || of free heavy quark decays. Next-to-leading order corrections, 
starting at 0(Xq CD /Mq), are expected to be small in heavy meson decays. This 
approach has also been extended to the case of nonleptonic decays |§ . 

To extract V u b from the B — > X u tv decay, one needs to measure the charged 
lepton spectrum near the high end of Eg, such that the huge B — > X c lv background 
stops contributing. The measurement must be performed within an accuracy of 
several hundred MeV, because the energy difference between the end points of the 
b — > c and b — > u transitions is only about 330 MeV. Unfortunately, OPE, the 
theoretical tool, breaks down in this region. Since the expansion parameter 1/Mq 
should be replaced by 1/(Mq — q • v) when q is not small at the end point, OPE 
is not reliable. The heavy quark velocity v is defined by Pq = Mqv, Pq being the 
heavy quark momentum. To circumvent the difficulty, Neubert 0] and Bigi et al. 



1 1 l]] have performed a resummation of OPE, that results in a model independent 
"shape function" in the description of the charged lepton spectrum. This universal 
shape function can be determined in principle by an infinite tower of nonperturbative 
matrix elements expanded in the increasing power of 1/Mq, and has been employed 
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in the study of the inclusive rare decay B — > X s 7 f[2 , I3[ . One can therefore measure 
this shape function, say, in the B — > X s 7 decay, and then apply it to the 5 — > 
decay to obtain a model independent prediction of V u b- However, under the present 
experimental situation, one has to make an ansatz for this shape function according 
to some QCD constraints. Data of the decay B — > will definitely remove this 
ambiguity in the choice of shape functions in the future. 



On the parton model side, Bareiss and Paschos [[L4J] argued that the phase space 
for the decays is dominated by distances near the light cone, and one can use the dis- 
tribution function of the Q quark in the heavy meson Hq in the infinite momentum 
frame to fill up the gap between the simple heavy quark kinematics and the heavy 
meson kinematics. This distribution function can be obtained by measuring the 
Q-quark fragmentation function in the heavy meson production from e + e~ annihila- 
tion. To remove the singularities near the end point of the charged lepton spectrum, 
which are due to soft gluon bremsstrahlung, these authors simply resummed soft 
gluon corrections in the naive leading (double) logarithmic approximation and ob- 
tained a Sudakov suppression factor. To fill up the kinematic gap, Altarelli and 



Petrarca |T3| regarded the light quark inside Hq as a quasi-free particle but with a 



Gaussian spectrum of Fermi-momentum p, 

/(p) = -|^exp(-^), (1) 

where pf is a free parameter that can be fixed by heavy quark symmetry. To 
smooth out the end-point singularities, these authors also resummed leading soft 
gluon contributions into a Sudakov form factor. 



Korchemsky and Sterman |16 gave the first PQCD treatment of the decays 
B — > X s 7 and B — > X u £u, in which higher-order corrections were factorized into a 
soft distribution function, a jet function and a hard scattering amplitude according to 
the kinematic regions of loop momenta. The equivalence between the shape function 
HlO| , 0] an d the distribution function was pointed out. Soft gluon corrections, 
which correspond to the nonperturbative origin of the distribution function, were 
resummed systematically up to next-to-leading logarithms using the Wilson-loop 
formalism. However, their analysis is appropriate only in the end-point region, and 
the effects of resummation were not estimated. 

In all the above approaches, the factorization was formulated at the quark-level 
kinematics, and the missing states inside the kinematic window between Mq and 
Mh, Mu being the Hq meson mass, were populated by introducing extra heavy 
quark Fermi motion, which arises from the recoil of the light partons in Hq. Our 
approach is formulated at the meson-level kinematics directly, in which the heavy 
quark Q, carrying a fraction of the heavy meson momentum, has an invariant mass 
close to Mq. Our formalism therefore removes the ambiguity in the definition of the 
heavy quark mass Mq, and the kinematic gap is filled up naturally. 

In this paper we shall derive the PQCD factorization formula for the semilep- 
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tonic decay B — ► X u £u in a rigorous way, which is suitable for the entire range 
of the spectrum. The factorization procedures demand the inclusion of the trans- 
verse degrees of freedom of partons. Hence, we perform the resummation of large 
perturbative corrections in the transversal configuration space using the technique 
developed in M, which is also accurate up to next-to-leading logarithms. It can be 



shown that our resummation result coincides with that in |16] at the end point. The 
transverse momenta carried by the b quark inside the B meson, whose distribution 
is governed by the Sudakov form factor from the resummation, play an important 
role here. They diminish the on-shell probability of the outgoing u quark at the end 
point, and thus suppress the singularities. 

We define the kinematics of the inclusive semileptonic decays of heavy mesons 
in Section 2, and derive the factorization formula for the charged lepton spectrum, 
which incorporates the transverse degrees of freedom of partons. The formula is 
expressed as the convolution of a hard scattering amplitude with a jet and a universal 
soft function. In Section 3, we resum the large logarithms in these convolution factors 
by solving a set of evolution equations. The initial condition of the soft function 
is identified as the distribution function, which is equivalent to the shape function 
mentioned above. In Section 4, we construct a distribution function according to the 
HQEFT based OPE, and relate the mean and width of the distribution function to 
the hadronic matrix elements of the kinematic operator. Hence, both perturbative 
higher-order corrections and nonperturbative corrections from Fermi motion are 
included in our formalism. We present numerical results in Section 5, and show 
that the Sudakov form factor from the resummation and the distribution function 
indeed render the end-point sprectrum smoother as stressed in [|12| . Section 6 is the 
conclusion. 



2 Factorization Theorems 

We consider the semileptonic inclusive decays of a B meson, 

B(P B ) -> £(jpt) + v(p v ) + hadrons . (2) 

The three independent kinematic variables are choosen as Ei, q 2 and go i n our 
discussion. E t and q = pe + p u have been defined in the Introduction, and q is the 
energy of the lepton pair. With these variables, the triple differential decay rate is 
written as 

d 3 T 1 

where Mb is the B meson mass, and the weak matrix element is given by 

M = v uh %l r> J2(x\r\B) . (4) 

V2 x 
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In eq. (£|) V u t, is the corresponding CKM matrix element, and = uT^b is the elec- 
troweak current with T M = 7 M (1 — 75). 



We work in the rest frame of the B meson, and choose the following light-cone 
components for relevant momenta, 

Pb = (Pb, p b , °±) , Pi = (Pe, °> 0-0 , Pv = (pt,P„, P»±) , ( 5 ) 

with Pg = Pg" = M b /a/2 and pi = 0. The independent variables are identified as 
p\ , p~ and p+, and their relations to E e , q 2 and g are P^ = p\ j q 2 = 1p\p~ 
and q = (pj + p+ + p~)/v2, respectively. We define p, = P# — p as the 6 quark 
momentum, which satisfies P 6 2 ~ M 2 , M& being the 6 quark mass, p is the kicks 
from the light components inside the B meson, which has a large plus component 
p + and small transverse components p_i_. The purpose of introducing the transverse 
degrees of freedom will become clear later. The b quark decays into a u quark with 
momentum P u = Pb — p — q. We have distinguished the B meson momentum Pg 
from the b quark momentum P& here. 

It is more convenient to employ the scaling variables 

2E e q 2 2q 

instead of the dimensionful ones Eg, q 2 and go- Note that the scaling variables are 
defined in terms of the B meson mass Mb, since we formulate the factorization 
according to the B meson kinematics. This differs from the conventional treatment 
in the literature where the scaling variables were defined in terms of the b quark 



mass. For massless leptons, it is easy to show, using the momentum configurations 
defined in eq. (|), that the phase space is given by 

0<x<l, 0<y<x, - + x<y <y + l. (7) 

x 

In the end point region with x — > 1 (pi — > Mb/V%) and y — > (p~ — > 0), we 
have yo — > 1 (p+ — > 0) and p — > 0. The u quark then has a large minus component 
P~ = (1— y /x)M B / V% but a very small plus component P+ = (l — y — y/x)M B /\/2, 
and thus a very small invariant P 2 = M B (1 — yo + y), which forms an on-shell jet 
subprocess. The u quark travels a long distance of C(1/Aq C d) before hadronized. 
Besides, the B meson is dominated by soft dynamics, which is the origin of the soft 
function stated in the Introduction. The remaining dominant subprocess is the hard 
one, which contains the weak decay vertices. Therefore, the important contributions 
are factorized into the soft (S), jet (J) and hard (H) subprocesses as shown in fig. 1. 

The factorization formula for the inclusive semileptonic decay B — > X u iv is 
written as 

- M% dz d 2 p ± 



H 0) dxdydy 
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xS{z,p±,fj l )J(z,P u ,p ± ,fj,)H(z,P u ,pj.,/x) 



with the momentum fraction z defined by z = P^~/P B = 1 — p + /Pg and vf^ = 
— —\V u b\ 2 Mg [ 16| . \i in eq. (§) is the renormalization and factorization scale. The 



167T 3 i 

triple differential decay rate is, of course, /i independent. Note that in the region 
y — > x ~ 1 the outgoing u quark becomes soft and eq. @ fails. We shall show that 
contributions from this dangerous region are suppressed by phase space. The upper 
limit of z takes the value z max = 1 in our analysis. If performing the factorization 
according to the b quark kinematics, one must assume z max = M B /M b , which is 
greater than 1, in order to fill up the kinematic window. It has been explained [[16] 



that z max > 1 is not allowed in perturbation theory, and is thus of nonperturbative 
origin. From the kinematic constraints in eq. (0) and the on-shell condition of the 
w-quark jet, the lower limit of z should be z m i n = x, instead of z min = 0. 

The tree-level expressions for the convolution factors J and H are given by 



J(°) = 6(P£) =5\M_ 



l-y + y-(l-z)(l--) — 2 —~2 



x M B M B 



H<® oc (P b -p„)(pe-P u ) = ((PB-p)-Pv)(pe-Pu 

y , 2p± jp„_ L 

t 2 B 



oc (x-y) \y -x-(l-z)- + — — 2 ). (9) 



Equation @ can be regarded as an expression at the intermediate stage in the 
derivation of conventional factorization theorems. If the p_i_ dependence in J and H 
is negligible, the variable p_i_ in S can be integrated over, and eq. (H) reduces to the 
conventional factorization formula. However, it is obvious from eq. @ that at least 
the p_i_ dependence in J is not negligible, especially in the end-point region. This is 
the reason we introduce the transverse degrees of freedom into our analysis. 

Suppose we consider higher-order corrections to eq. @ from a gluon crossing 
the final state cut, and route the loop momentum £ through, say, the jet subpro- 
cess. Without losing generality, we approximate J{p + + £ + ,Pu + £~ ,p± + i±) ~ 
J(p + ,P~, p_i_ + £±) according to the kinimatic relations £ + < p + , £" < P~ and 
£±_ ~ P±- Hence, the loop integral cannot be performed unless the dependence of 
J on transverse momentum is known. This difficulty can be removed by Fourier 
transform, 

J(P + , Pu, P± +*±)=J tS* j(P + , Pu, h)e^-> h , (10) 

where the impact parameter h (Fourier conjugate variable of p±) measures the trans- 
verse distance travelled by the jet. Using eq. (|10|), the £± dependence is decoupled 
from the jet function, and the factor e t£± ' h is absorbed into the loop integral, which 
can then be performed. Therefore, an extra factor e l£j -' h is associated with each 
gluon crossing the final state cut in our formalism. 
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To further simplify the factorization formula, we neglect those terms involving 
p„± in J and H. This is a good approximation for x — > and x — > 1, since for 
x — > contributions from transverse momenta are not important, and for x — > 1 the 
magnitude = Mb\J (yo — x — y/x)y/x vanishes. Equation (§) then becomes 

1 d 3 T 2 yi r d 2 b 

Tfdxdydy ~ B L 

xS(z,b,ii)J(z,p-,b,fi)H(z,p-,fi) . (11) 



3 Resumming the Jet, Soft and Hard Subpro- 

cesses 

It can be shown that the dominant subprocesses J, S and H contain large loga- 
rithms form radiative corrections. In particular, J gives rise to double (leading) 
logarithms in the end point region. These large corrections spoil the perturbation 
theory and must be organized. In what follows we shall demonstrate in details how 
to resum these large corrections up to next-to-leading logarithms. The first step in 
resummation is to map out the leading regions of radiative corrections. We work 
in axial gauge n ■ A = 0, where n is a space-like vector. The 0(a s ) diagrams that 
contain large double logarithms in axial gauge at the end point are figs. 2a and 2b. 
The self-energy diagram fig. 2c and the diagram with a soft gluon connecting the 
two heavy quark lines, as shown in fig. 2d, give only single soft logarithms. 

In the collinear region with the loop momentum i parallel to P u and in the soft 
region with i — > we can eikonalize the heavy b-quark line. Then figs. 2a and 2b 
are factorized out of the cross section, and they are the diagrams that J absorbs, as 
shown in fig. 3. With the eikonal approximation, the b-quark propagator is expressed 
as 1/v • i to the order 1/Mb with v = (1, 1, 0j_). Hence, the factorization in eq. (]TTf) 
is in fact valid to 0(1/ Mb)- The physics involved in this approximation is that a 
soft gluon or a gluon moving parallel to P u can not explore the details of the b quark, 
and its dynamics can be factorized. This is consistent with the HQEFT, where the 
b quark is treated as a classical relativistic particle carrying color source. Since the 
large mass Mb does not appear in the eikonal propragator, the only large scale in 
J is P-. 

The remaining diagrams figs. 2c and 2d that give single soft logarithms are 
grouped into S. It is then obvious that S depends only on the properties of the 
bound state \B), but not on the particular short distance subprocess. Therefore, 
S is a universal function describing the distribution of the b quark inside the B 
meson. In fact, in the 1/M& limit one can identify S as the model independent 
"shape function" or "primodial function" obtained from the resummation in [0, |13] 
and |jTI|| , respectively. 
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The basic idea of the resummation technique is as follows. If the double log- 
arithms are organized into an exponential form, J ~ exp[— lnP^ ln(lnP~/ln6)], 
one can simplify the analysis by studying the derivative of J, dJ/d\nP~ = C J 0. 
Since the coefficient function C contains only single logarithms, it can be treated 
by renormalization group (RG) methods. In this way, one reduces the complicated 
double-logarithm problem to a single-logarithm problem. 

J is scale invariant in the gauge vector n as shown by the gluon propagator in 
axial gauge, -^r^p , with 



N^(£) = gT - , + n 2 -— - 2 . (12) 

n ■ t [n ■ ty 

Therefore, J must depend on P~ through the ratio (P u ■ n) 2 /n 2 . It is then easy to 
show that the differential operator d/d In P~ can be replaced by d/dn using a chain 
rule, 

-f- = -^-v'° J-J (13) 
d In P~ v 1 ■ n dn a 

with the vector v' = (0, l,0j_). This simplifies the task tremendously, because the 
momentum P u flows through both quark and gluon lines, but n appears only in 
gluon propagators. 

Applying d/dn to the gluon propagator, we obtain 

d -N^ = _J_(jVT + N va t) . (14) 



dn a t ■ n 

The momentum I that appears at both ends of the differentiated gluon line is con- 
tracted with the vertex, where the gluon attaches the u quark or the eikonal lines. 
Next we add up all diagrams with different differentiated gluon propagators, and 
apply the Ward identity. Finally, we arrive at a differential equation for J as shown 
in fig. 4, where the square vertex represents 

n 2 

gT a ( , M9 / ° , (15) 

[v ■ n){i ■ n) 

with T a the color matrix. The factor 2 counts the two external quark lines of J. An 
important feature of the square vertex is that the gluon momentum i does not give 
rise to collinear divergences because of the nonvanishing n 2 . The leading regions of I 
are then soft and ultraviolet, in which the subdiagram containing the square vertex 
can be factorized as shown in fig. 5 at 0(a s ). Hence, the differentiation really turns 
the double-logarithm problem into a single-logarithm problem as stated before. 

To seperate the soft and ultraviolet scales in fig. 5, we introduce a function K, 
to organize the soft logarithms from the four diagrams figs. 5a-5d, and a function Q 
for the ultraviolet divergences from figs. 5e and 5f. Double counting is avoided by 
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the subtraction in Q. Generalizing /C and Q to all orders, we derive the differential 
equation for J, 



d 



d\nP: 



-J(P-,M= 2 K,{bn,a a (n)) + G(P^ln,a a {n)) J(P-,M- (16) 



The scale b in /C arises from fig. 5b, which contains an extra factor e t£± ' h as explained 
in Section 2. 

At the one-loop level, the renormalized function /C calculated in dimensional 
regularization for the gauge vector n oc (1, —1, 0j_) is given by, 



K. = fig. 5a + fig. 5b + fig. 5c + fig. 5d - 5K , 



with 



fig. 5a + fig. 5b 



(17) 



-9 2 C F ^ J 



d 



4-e/ 



(2tt) 4 - 

'2 



n 2 v' a v' p 



a s 
2^ 



£ 2 + ie) (v' ■ n)(£ ■ n){£ ■ v') 
Cf^ + Invr /i 2 6 2 e 7 | , 



-N^{£) 



(18) 



fig. 5c + fig. 5d = -g 2 C F fi e J 



(2k) 

= . (19) 

8K, is the scheme- dependent counterterm which will be specified later. C F = 4/3 is 
the color factor, and 7 is the Euler constant. Similarly, Q is given by 

G = fig. 5e + fig. 5f - 5Q , (20) 

with 

d*-*£ p - a 



fig. 5e + fig. 5f = -g 2 C F fi £ J 

r p. 



(27r) 4 " e (P u -£) 2 + ie 



2„,/ 



n v 



(3 



N<*P(£) 



£-v' a ) (v'-n)(£-n) £ 2 + ie 

lH I e TTjd z 



-^1-: + "^^ ■ Pi) 



and z/ = (n • v') 2 /\n 2 \ is the gauge factor. 
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In the MS scheme with 

SJC = sg = - 

the one-loop /C and Q are given by 



^C F (-+ In 4^ 
Z7r V e 



(22) 



/c = 




(b 2 ji 2 e 2 ^ 






V 4 


= 










V /^ 2e 



(23) 



We have seperated the two scales in J, 6 and P~, into the functions K and Q, respec- 
tively, such that RG methods are applicable to the summation of the corresponding 
single logarithms. Although K, and Q possess individual ultraviolet pole, their sum 
is finite, and thus a RG invariant quantity The RG equations for the renormalized 
K and Q are 



with 



H— /C(6/x, a s {n)) = -7K(a«(^)) = -^Q{ P u /v,<Xs(p)) 
djj djj 



Ik = i^-^-SJC = -ji-^-SQ 

djj djj, 



the anomalous dimension of /C. To two loops, 7^ is given by [|Tj 



a 



Ik = —C F + — C 



71 



7T 



/67 vr 2 \ 5 

64 n, 

1 36 12/ 18 3 



(24) 



(25) 



(26) 



where rif = 4 is the number of quark flavors and Ca = 3 is the color factor. 
Equation (|24D has the solution 

IC(bji,a s (ji)) + g(p-/ji,a s (jj)) = K(l,a s (p-)) + g(l,a s (p-)) 

p «~ dp . . , . . 
— A(a s (/i)) 

1/6 /Li 



with the anomalous dimension 



d 



A(a s ) = 7^(a s ) + (3(g)— a s ) . 

Substituting eq. (^7|) into Ql6]), we obtain the evolution of J in P~ and 6, 

J(p-,b,jj)=exp[-2s(p-,b)]J(b,jj) . 



(27) 



(28) 



(29) 
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The RG invariant Sudakov exponent is given by [JT8 



gin I) + 



A (2) 



2A ' V /, 



4ft 2 U 



1-1 



^4(i) 
2A 



(g - S) 



-g 



ln(26) + 1 ln(2g) + 1 



4/3? 

^(2) A (D 



4/3 2 4A 



In 



,27-1' 



In 



with the variables 



' 8/3? 



g = In 



ln 2 (2g) - ln 2 (26) 



P. 



u 

A" 



In 



bA 



(30) 



(31) 



The QCD scale A = Aqcd will be set to 0.2 GeV in the numerical study in Section 
5. The coefficients /3j and are 



01 

02 

A® 

A (2) 



33 -In 



f 



12 ' 
153 - 19n/ 

24 

4 

3 ' 

67_ 7T> 
9 3 



10 8 ^ , 
— Wf H — pi In 

27 ; 3 M 



(32) 



Having summed up the double logarithms, we concentrate on the single loga- 
rithms in S, H and the initial condition J(b,fi). Since both the differential decay 
rate and the Sudakov exponent s(P~, b) are RG invariant, we have the following RG 
equations: 

DJ(M = -2 lq j{b,(i), 
T>S(b,fi) = -7 S 5(6,/i), 
VR(p-,ii) = (2 lq + ls )H(P;,n), (33) 

with 

D = „| + (34) 

7 g = —a s /iT is the quark anomalous dimension in axial gauge, and 75 = — {a s /n)CF 
is the anomalous dimension of S. The function S in fact contains soft single loga- 
rithms only from fig. 2d, because the contribution from fig. 2c vanishes for the guage 
vector n oc (1, — 1, 0j_) under eikonal approximation. Hence, 7s is derived from the 
evaluation of fig. 2d straightforwardly. 
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Integrating eq. 
3(z,b,fi) 



we obtain the evolution of all the convolution factors, 



exp 



exp 



exp 



f 1 dp, 



, —IfqMn)) 
b /i 



J(z,b,l/b) 



1/6 (1 



f(z,b,l/b) 



— [27,(a fl (/x)) + 7s(a s (/i))] 
/i 



H(z,P-,P-).(35) 



We shall neglect the intrinsic b dependence of the distribution function / below. If 
Sudakov suppression in the large-6 region is strong, we may drop the evolution of / 
and J in b, which is proportional to a s (l/b). Hence, we assume f(z,b, 1/b) = f(z), 
J(z, b,l/b) = j(°'(z, b), the Fourier transform of the tree-level expression in eq. (|9]), 
and H(z,P-,P-) = H(°\z,P-). 

Substituting eq. (^) into (|TT|), we derive the factorization formula for the inclu- 
sive semileptonic B meson decay, 

1 d 3 T 



rj 0) dxdydy 



Mlfdz r b ^f(z)J^(z,b)H i0) (z,P- 
Jx Jo 2n 



xexp[-S(p-,b)}. 
The complete Sudakov exponent S is given by 

S(P~,b)=2s(p-,b) 



m b 



(36) 



(37) 



with P~ = \n(P~/A), which combines all the exponents in eq. and includes both 
leading and next-to-leading logarithms. It is straightforward to observe from eq. (|37j) 
that the Sudakov form factor e~ s falls off quickly at large b ~ 1/A, where a s {l/b) > 1 
and perturbation theory fails. Hence, the Sudakov form factor guarantees that main 
contributions to the factorization formula come from the small b, or short- distance, 
region, and the perturbative treatment is indeed self-consistent. 

We stress that our formalism is applicable to the entire range of the spectrum, if 
the dependence in J and H was not neglected, not only to the end-point region 
as in JT^]. We neglected the p u ± dependence for the sake of simplicity of calculation. 
Since the Sudakov form factor is defined in the region 

M B „_ 1 



1 



V2 



p: 



>->A 



(38) 



there exist the suppression effects as y < (1 — \^2A/Mb)x. The phase space with 
suppression expands to the largest extent for x — > 1 and vanishes at x — ► 0. Hence, 
lowest-order predictions receive maximal suppression at the end point, correspond- 
ing to the presence of large logarithms, but are modified only slightly at small x, 
implying weaker logarithmic corrections. This is consistent with our expectation for 
the Sudakov effects in the whole range < x < 1. 
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4 Constructing the Universal Soft Function 



Since all the double logarithmic corrections have been absorbed into the jet 
subprocess, the soft function S contains only soft single logarithms from fig. 2d 
(fig. 2c does not contribute because of the choice of the guage vector n oc (1, —1, OjJ). 
These single logarithms can be summed by solving the RG equation VS = —jsS as 
shown in eq. The solution has been given in eq. Q55|), 



S(z, b, /i) = exp 



, -7s « s (w) 

1/6 \l 



(39) 



where the initial condition f(z) for the evolution of S must be determined phe- 
nomenologically. From the definition of S, it is obvious that / depends only on 
the properties of the bound state \B), but not on the particular short distance sub- 
process. Therefore, / is a process-independent universal function describing the 
distribution of the b quark inside a B meson. 

As mentioned in the previous section, in the 1/M& limit one can identify / as 
the model independent shape function or primodial function obtained from the re- 
summation of OPE in Refs. |10| and [|ll| [T^j. It was observed by Korchemsky 
and Sterman |TB| that one can deduce some nonperturbative information of / from 
the perturbative resummation of soft gluons. Theses authors argued that the in- 
frared renormalons appearing during the resummation procedure produce power- 
correction ambiguities. Hence, there should exist corresponding nonperturbative 
power-correction ambiguities to render the physical / well defined. To leading or- 
der, with a "minimal" ansatz for integrating the first infrared renormalon singularity, 
Korchemsky and Sterman ITBfl derived 



fK S (z) 



V2 



=exp 



(1-*)' 
2a 2 



(40) 



which describes a Gaussian distribution with width a around z = 1. Requiring that 
the minimal ansatz be consistent with HQEFT and OPE, one fixes the width to be 



a 



j4_ 

3M 2 ' 



(41) 



where /i 2 = 0.54 ± 0.12 GeV 2 |L9| is obtained from QCD sum rule estimation. 



Recall that fxs was derived by integrating only the first infrared renormalon. 
This fact is reflected in the nonvanishing property of fxs m the unphysical region 
z > 1, and in its validity in the vicinity of z — 1. However, comparing fxs with the 
phenomenological distribution function [OT, 



fp(z) = N P 



z(l-z) 



[(1 



zY + e P z] 2 ' 



(42) 
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one finds that they do share the same leading (1 — z) 2 behavior, fp was obtained 
from the experiments of the B meson production in e + e~ annihilation by applying 
a crossing to the light quark and then a time reversal transformation p0| . The 



constant Np = 0.133068 is the normalization, and ep = 0.006 is the shape parameter. 
Inspired by this observation, we postulate the following two-parameter distribution 
function for the B meson, 

^ V-VV ' 1 -*'- (43) 

which has the correct leading (1 — z) 2 behavior near z — 1. The purpose of intro- 



ducing one more parameter a in eq. ( 4"3"|) is to allow a consistency check on fp in 
eq. (H). 

We now relate the parameters a and e to the hadronic matrix elements of some 
local operators derived from HQEFT based OPE, which are standard techniques 



TO, 12, 13]. The distribution function in axial gauge is defined by 



/*(*) = ^ / ^T el(1_2)P> ~ (B\K(0)b v (y-)\B) , (44) 

in which the large momentum of the rescaled heavy quark field b v (x) has been 
projected out as usual in HQEFT by 

b v (x) = e iMbV - x b(x) . (45) 

Hence, is independent of the 6-quark mass as it is written in terms of a 

matrix element in eq. (0). Note that the Dirac matrix structure has been properly 
factorized as shown in eq. (H). 



To connect with HQEFT, we write the b quark momentum as P^ = Pp — p = 
Mf)V + (Mp — Mb)v — p, and identify the residual momentum of the b quark as 
k = Av — p, A = Mp — Mb being the effective mass of the light partons in the B 
meson. The probability to find a b quark with light-cone residual momentum k + 
inside the B meson, f r {k + ), has been defined in |lj|. The moments of f r (k + ), 



A n = J k +n f r (k + )dk + } (46) 
have been derived in HQEFT based OPE as [O 

M 2 

Ao = l, A 1 = 0, A 2 = -^-K b ,... (47) 

These moments are expressed in terms of hadronic matrix elements corresponding 
to the structure of the 1/M& expansion, with 
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The vanishing 0(Aqqd/Mb) contribution to the first moment A\ is consistent with 
the conclusion from the renormalon analysis that the first nontrivial power correc- 
tion begins at 0(A 2 QCD /M B ) and with our intuition for vanishing; average residual 
momentum of the b quark inside the B meson in the heavy quark limit. 

The moments of /b(z) can be expressed as local hadronic matrix elements by 
performing an OPE of the bilocal operator b v (0)b v (y~) in eq. in the power of 
1/Mfe. The relation between fg and f r is then given by |13 



f B {z)dz = [f r (k + ) + 0(1/M b )\ dk + , (49) 
which reflects the difference of order A between the B meson kinematics and the b 



quark kinematics. Using eq. fl46|) and the definition z = P t ,/Pb or k + = Pb(z 



Mb/ Mb), it is straightforward to derive the moments of f B - They are 

[ 1 f B (z)dz = l, (50) 
J o 

Jl dz(l - z)f B (z) = A/ Mb + 0(A 2 QCD /M 2 B ) , (51) 

jf dz(l - zff B {z) = ^T + l K b + 0(A 3 QCD /M B ) . (52) 

The first formula gives the correct normalization of f B , which corresponds to the 
total number of b quarks inside a B meson. The second formula is related to the 
effective mass of light quarks, A. The third formula gives the hadronic matrix 
element of the kinematic operator, K b . 

To have a better insight, we examine if the distribution function is consistent 
with our physical intuition for the behavior of the heavy b quark inside a B meson. 
We calculate the mean fi and the variance a 2 of /b(^) from eqs. (|5"ID and (0), and 
derive 



^-w B + °(^k ] > (53) 

^ + c /a|A (54) 



Substituting the QCD sum rule [[L9j and B* — B mass splitting 0] results, 



M B = 5.279GeV , M b = 4.776GeV , K b = 0.012 ± 0.0026 , (55) 

we obtain fi = 0.90 and a 2 = 0.0080 ± 0.0017, implying that /b(z) peaks sharply 
around z«/i«l and has a width of 0(Aqct>/Mb). The parameters N, a and e in 
Jb(z) can be determined using eqs (|57j|)-(|5"2"|) with (|55| ) inserted, which are 

N = 0.02609, a = 0.9752, e = 0.001699 . (56) 

The value of a, derived from the QCD constraints (by taking finite number of mo- 
ments only), is very close to unity. This is consistent with the expectation from fp. 
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However, we emphasize that fp is not quite consistent with HQEFT, its first and 
second moment differing from eqs. ( [51] ) and ( [52] ) by at least 45%. 



The distribution functions in eqs. (H2T) and (E3|) will serve as the initial conditions 



of the soft function in fl39|). We then derive the PQCD factorization formula for the 
inclusive decay B — > X u £u, with the phenomenological inputs satisfying the QCD 
constraints from HQEFT based OPE. 



5 The Charged Lepton Spectrum 



In this section we evaluate eq. fl36|) numerically for various distribution functions. 
The charged lepton spectrum for the decay B — > X u lv from the naive quark model is 
obtained by simply choosing f(z) = 5(1 — z) and ignoring the transverse momentum 
dependence in and H^°\ A simple calculation leads to 

1 dT x 2 , 

(3 - 2s) , (57) 



rV 0) dx 6 

which corresponds to the solid curve (1) in fig. 6. This curve does not fall off at 
the end point of the spectrum, contradicting the observed behavior of the inclusive 
semileptonic decays of B mesons. The discrepancy implies that the tree-level anal- 
ysis is not appropriate, especially in the end-point region where PQCD corrections 
are important as discussed in Section 3. 

We then take into account Sudakov suppression from the resummation of large 
radiative corrections. Substituting f(z) = 5(1 — z), = (x — y)(yo — s) and the 
Fourier transform of J (0) = 5(P 2 ) with P 2 = M|(l -y + y-p 2 ± /M^) into eq. fl3J), 
we derive the modified quark model spectrum. This spectrum is, after integrating 
eq. (B9) over z and yo, described by 



1+y- x)J l (r ] M B b) - -^-r)J 2 (vM B b) + V 2 J 3 ( V M B b) 

M B o 

where P~ = (1 — y/x)M B /y/2, r\ = (x — y)(l/x — 1) and Ji,J 2 ,J 3 are the Bessel 
functions of order 1, 2 and 3, respectively. Note the presence of the Sudakov form 
factor e~ s and the expression in the square brackets which comes from J(°\ The 
cutoff 1/A of the impact parameter b is set by the Sudakov form factor. Numerical 
results of eq. (|58|) for A = 0.2 GeV are shown by the dashed curve (1) in fig. 6. 
Since we have neglected the p u ± dependence in J and H for simplicity, eq. (|58| ) 
is appropriate only for small and large x. Therefore, to obtain the dashed curve 
(1), we evaluate eq. fl58|) in the regions < x < 0.7 and 0.9 < s < 1, and then 
extrapolate from x = 0.7 to 0.9 smoothly. The dependence on A in our analysis is 
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also examined, and it is found that predictions increase by only 10-20% if A was set 
to 0.1 GeV. 



One observes immediately that the Sudakov effects alone are enough to render 
the uprising free quark spectrum fall off at the end point. This is consistent with 
our expectation that the inclusion of transverse momenta and Sudakov suppression 
diminishes the on-shell configuration of the outgoing w-quark jet. Another important 
feature in fig. 6 is that the solid and dashed curves coincide with each other in the 
region x — > 0. This indicates that the Sudakov effects almost cease to contribute 
away from the end point as stated in Section 3. 

The spectrum from the parton model without Sudakov suppression is obtained by 
adopting H<® = (x-y)[y -x-(l-z)y/x] and P u 2 = M%[l-y +y-(l-z)(l-y/x)]. 
With integration over y , we derive 

1 dT f x r l 

dy dzf(z)(x-y)(y + z-x) , (59) 



r/°> dx 



o 



where f(z) can be replaced by the distribution functions given in eqs. fl4"2"|) and 
(f43"l). Predictions from the use of fp and fp are represented by the solid curves (2) 
and (3) in fig. 6, respectively. Both of the spectra deviate from the quark model 
one slightly at small x, and vanish at the end point. Since eq. ( |59| ) incorporates 
nonperturbative effects from primodial heavy quark motion [O, [T3j (or soft dynamics 
in our formalism) through /, we conclude that these nonperturbative corrections are 
indeed important in the end-point region. 

At last, including Sudakov suppression into eq. (H), we arrive at the charged 
lepton spectrum of the B — > X u iv decay that takes into account both large pertur- 
bative and nonperturbative corrections, 



1 dT r x /"VA r 1 



M B b 



2 iHiMBb)+emMBb) 



{z + y-x^ZMBb) 
s(Pu,b) ; (60) 



with £ = J (x — y)(z/x — 1). Predictions from fp and fp are shown by the dashed 
curves (2) and (3), respectively. They coincide with the solid curves at small x, but 
descend by about 50% at x — > 1 as shown in fig. 7, implying strong suppression in the 
end-point region. The slope of the spectrum then becomes smoother as expected. 

From fig. 6 we evaluate the total decay rate T/T^, and results along with the 
Sudakove effects are displayed in Table. 1. We find that the overall suppression 
from the Sudakov effects is 8% for the quark model and less than 7% for the use 
of fp and fp. The two distribution functions lead to about 20% difference in the 
total decay width. Comparing to the drastic distinction between fp and fp with 
Np/N ~ 5 and ep/e ~ 4, our formalism is quite insensitive to the choice of different 
distribution functions. The overall Sudakov suppression of less than 8% indicates 
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that PQCD corrections are actually not important for most part of the spectrum. 
This is consistent with the fact that corrections from transverse momenta are an 
0(1/M 2 Q ) effect. 

Note the 30% suppression on the quark-model results from the distribution func- 
tion fp (the suppression from fp is even stronger). It has been found in HQEFT 
that effects from nonperturbative corrections are only of 0(1/ M B ), which should be 
less than 5% (see p2] and references therein). The small nonperturbative correc- 



tions to the total decay rate are closely related to the vanishing first moment of the 
residual momentum structure function f r . This apparent discrepancy can be traced 
back to the fact that the B meson kinematics is employed in our formalism, while 
the 6-quark kinematics is employed in the conventional approaches. Hence, in our 
quark-model analysis the b quark in fact carries the full momentum p, = Pp, and 
thus the charged lepton energy Ei can reach the maximum Mp/2 (x = 1). This 
momentum configuration is allowed in factorization theorems if transverse degrees 
of freedom of partons were included, because its invariant P 2 = M\ — may still 
be close to the mass shell in the region without Sudakov suppression. Without kp, 
eq. ( [57] ) should be regarded as an expression that is generated in our formalism to 
bear the same form as the leading-power results in HQEFT. For a free b quark with 
momentum M^v, Eg can only reach M&/2, instead of Mp/2. Strickly speaking, our 
quark-model predictions and the leading-power predictions in HQEFT have different 
meanings. 

We stress that our results do not violate the conclusion from HQEFT, if they 
were interpreted in a proper way. To confirm this, we identify p, = (M 2 /2P B , P B , 0) 
as the momentum carried by a free b quark in factorization theorems for B meson 
decays, where the minus component P b ~ has been set to P B . That is, the free b 



quark is not at rest inside the B meson. We then reexpress eq. (|59 ) into a form 
similar to that in 



1 F(x) o(MI- x \ + f(^\m( x) , (61) 



T^dx v ' \M 2 ) \M 2 B/ 

with F(x) = x 2 (3 — 2x)/6 being the quark- model prediction derived from the con- 
ventional approaches and 

F l-^) M(x) = j o dy j x dzf(z)(x -y)( y + z -x)- F(x)6 j-| - . (62) 



The step function in eq. (pi ) specifies the maximal Eg in the decay of a free b quark 
with the above momentum P&. The function M(x), representing nonperturbative 
corrections to the b quark decay, coincides with the shape function S(x) defined in 



22 . 



We shall show that the contribution from M(x) to the total decay rate is indeed 
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of 0(1/M B ). Integrating eq. fl62| ) over x, we obtain 

F \W B ) I M{x)dX =12/0 dZz4f{z) ~ Jo • (63) 



An arbitrary structure function /, which possesses the same moment as in eq. (51), 
can be expanded in terms of ^-functions: 

f(z) = 6(1 -z)- A^i - Z ) + 0(A 2 /Ml) . (64) 



Inserting eq. (|64]) into (|63|), we justify straightforwardly that the nonperturbative 
correction 



1 A ' A 0(A 2 /M 2 B ) (65) 



3 M B 3 M b 



vanishes at 0(1/ Mb) as concluded in f22fl . In summary, the quark-model contribu- 
tion from the window between x = M^/M B and x — 1, with a width of 0(1/ Mb), 
cancels the 0(1/ Mb) correction from the structure function, such that the nonper- 
turbative correction is of 0(1/ M B ). 



According to eq. (|63D, the suppression from nonperturbative corrections is about 
5% for fs and more than 10% for fp. The percentage for fp is still large, because its 
first moment does not satisfy the requirement of HQEFT, and thus the cancellation 
at the power 1/M B is not complete. 

The ambiguity from the choice of distribution functions can be removed, once 
the spectrum of the decay B — > X s 7 is available. We can fix the universal B meson 
distribution function from these data, substitute the distribution function into our 
formula, and predict the end-point spectrum of the decay B — > X u tv. A model 
independent extraction of the CKM matrix element \V u b\ then becomes possible 



6 Conclusion 

We have studied the inclusive semileptonic B — > X u tv decay using the PQCD 
formalism. In order to seperate the B X u £u signals from the B — > XJLv back- 
ground, we must investigate the charged lepton spectrum near the end-point region 
within an accuracy of about 330 MeV. It has been found that there exist large 
perturbative corrections in this region, which are resummed into the Sudakov form 
factor and included into the factorization formula. The transverse degrees of free- 
dom of the b quark diminish the on-shell configuration of the outgoing u-quark jet. 
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The quark- model spectrum then falls off at the end point, consistent with the ex- 
perimental observation. There is no ambiguity associated with the kinematic gap, 
because we formulate the factorization for B meson, instead of b quark, decays. 

We have constructed a distribution function, whose parameters are determined 
by the HQEFT based OPE, and whose width and mean are related to hadronic 
matrix elements of the kinematic operator. These hadronic matrix elements are then 
fixed by QCD sum rule results ]19| and B* — B splitting data The distribution 



function, absorbing important nonperturbative corrections from heavy quark Fermi 
motion, can also render the quark model spectrum vanish at the end point. 

We emphasize that our formalism incorporates both large perturbative and non- 
perturbative corrections in the end-point region of inclusive B meson decays in a 
systematic way, and that it provides a natural normalization for the spectra. This 
enables the direct extraction of the CKM matrix element \V u b\ from experimental 
data. When more data are available, our formalisms can also be used to test PQCD 
in B meson decays. 



It is an important issue that current experimental data [23| of the B — > X£v 
branching ratio suggest 

BR(B -> Xlv) < 11% . (66) 

The naive quark-model prediction for this branching ratio is more than 15%. Al- 
though PQCD suppression at the end point is around 50%, the overall suppression 
amounts to 8% at most. With modification from the massiveness of the charm quark, 
our formalism can be applied equally well to the semileptonic decay B — > X c iv. 
Hence, we conclude that PQCD corrections suppress the overestimated theoreti- 
cal value of the semileptonic branching ratio only down to 13.8% at best. On the 
other hand, the distribution functions may decrease the quark-model predictions 
by about 30%. However, the distribution function is universal as stressed before, 
and it is very plausible that it gives an equal amount of suppression to nonleptonic 
decays. Therefore, introducing a distribution function may not be able to remove 
the disagreement. 

Based on the above discussion, we propose three possibilities to resolve the dis- 
crepancy: (1) the distribution function suppresses semileptonic B meson decays 
maximally, but does nonleptonic decays minimally. This may arise from the dif- 
ferent phase space in these two cases. (2) Factorization theorems break down in 
B meson decays. (3) New QCD effects or new physics appears. Blok and Mannel 



24|] argued that factorization theorems may still hold, and thus the confrontation 
between data and theoretical predictions becomes acute. To settle down the issue, a 
careful PQCD analysis of B meson nonleptonic decays is required. We shall discuss 
these subjects in a forthcoming article. 
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Table. 1 The total decay rates for the quark model and for the use of the distri- 
bution functions fs and fp. 



r/rf 


8(1 -x) 


!b 


fp 


without suppression 


0.0833 


0.0586 


0.0446 


with suppression 


0.0767 


0.0548 


0.0425 


Sudakov effects 


7.92% 


6.48% 


4.71% 
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Figure Captions 



Fig.l: 

Factorization of inclusive semileptonic decays of the B meson into a soft (S), a jet 
(J) and a hard (H) subprocess. 



Fig.2: 

Lowest-order radiative corrections to the inclusive B meson decays. 



Fig.3: 

Factorization of the jet subprocess. 



Fig.4: 

Graphic representation of eq. flilf ). 



Fig.5: 

Lowest-order diagrams for the functions /C and Q. 



Fig.6: 

Charged lepton spectra of the B — > X u tv decay for (1) f(z) = 5(1 — z), (2) f(z) = 
fB^z) and (3) f(z) = fp(z). The solid (dashed) curves are derived without (with) 
Sudakov suppression. 



Fig.7: 

Charged lepton spectra of the B — > X u lv decay near the end point for the use of 
fs and fp. Conventions are the same as those in fig. 6 but with the dotted curve 
corresponding to the dashed curve (3). 
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